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ABSTRACT 

Very metal-poor stars are of obvious importance for many problems in chem- 
ical evolution, star formation, and galaxy evolution. Finding complete samples 
of such stars which are also bright enough to allow high-precision individual 
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analyses is of considerable interest. We demonstrate here that stars with iron 
abundances [Fe/H] < —2 dex, and down to below —4 dex, can be efficiently 
identified within the Radial Velocity Experiment (RAVE) survey of bright stars, 
without requiring additional confirmatory observations. We determine a cali- 
bration of the equivalent width of the Calcium triplet lines measured from the 
RAVE spectra onto true [Fe/H], using high spectral resolution data for a sub- 
set of the stars. These RAVE iron abundances are accurate enough to obviate 
the need for confirmatory higher-resolution spectroscopy. Our initial study has 
identified 631 stars with [Fe/H] < —2, from a RAVE database containing ap- 
proximately 200,000 stars. This RAVE-based sample is complete for stars with 
[Fe/H] < —2.5, allowing statistical sample analysis. We identify three stars with 
[Fe/H] < —4. Of these, one was already known to be 'ultra metal-poor', one is 
a known carbon-enhanced metal-poor star, but we obtain [Fe/H] =— 4.0, rather 
than the published [Fe/H]=-3.3, and derive [C/Fe] = +0.9, and [N/Fe] = +3.2, 
and the third is at the limit of our S/N. RAVE observations are on-going and 
should prove to be a rich source of bright, easily studied, very metal-poor stars. 

Subject headings: stars: abundances — stars: Population II — Galaxy: abun- 
dances — Galaxy: stellar content 



1. Introduction 



Very metal-poor (VMP) stars, defined conventionally as having [Fe/H] < —2 ( iBeers &: Christlieb 



20051 ). are of interest in many outstanding problems in star formation and galaxy evolution. 
For example, the detailed shape of the metal-poor tail of the metallicity distribution of field 
stars of the Milky Way can distinguish different models of chemical evolution, constraining 
gas flows, the nature of t he progenitor systems in which stars form, and pregalactic enrich- 
ment (e.g. |Prantzosll2009l ). Elemental abundance patterns of VMP stars are sensitive to the 
high-redshift stellar Initial Mass Function, star-for mation process, and poss ibly identify the 
sources responsible for reionization of the Universe ( iBromm &: Larson! 120041 ) . 



Detailed elemental abundances of metal-poor stars require high-resolution high-quality 
spectroscopy. Thus apparently bright targets with high probability of being truly very metal- 
poor are desirable. We demonstrate in this paper that the RAVE survey is an ideal source 



1 Based in part on observations collected at the European Organization for Astronomical Research in the 
Southern Hmcisphere, Chile, in the framework of proposals 081.B-0900 and 080.B-0927 
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of bright metal-poor stars; these can be identified directly from the RAVE data, with no 
need for confirmatory higher-resolution spectroscopy. We derive the calibration by which 
the RAVE parameters may be placed on a true iron abundance scale, valid for [Fe/H] values 
as metal-poor as —4 dex, into the regime of 'ultra metal-poor stars' in the nomenclature of 
Beers & Christlieb (2005). 



2. Identifying Very Metal-Poor Stars in the RAVE Survey 
2.1. The RAVE Survey 



The Radial Velocity Experiment (RAVE) f lSteinmetz et al.ll2006l : IZwitter et al.ll2008f ) is 
a spectroscopic survey of apparently bright stars, / < 13, with some 350,000 stars observed 
to date. RAVE uses the six-degree-field (6dF) multi-object, fibre-fed spectrograph on the 
UK Schmidt telescope. Spectra of ~ 100 targets are obtained simultaneously. The dominant 
target population belongs to the thin disk, but includes many thick disk and halo giant 
stars. The spectra are centered on the infrared Calcium triplet, cover the wavelength range 
region AA8410 — 8795A and have resolving power R ~ 7500. This is significantly higher 
than that which is typically employed in surveys to fi nd metal-poor candida tes for later 



confirmation - for example, the Hamburg/ESO survey (IChristlieb et al.ll2004l ) is based on 
objective-prism spectra with R ~ 400 at the Ca II K line at 3934A. We show here that 
RAVE has the resolution to determine the metal-poor nature of a star directly from the 
survey spectra, without the need for follow-up confirmation, providing ideal targets for later 
detailed analysis. 

The RAVE spectral analysis pipeline provides estimates for each star of the line-of- 
sight velocity and the stellar parameters (T e g, log g, [m/H], and v sini) through \ 2 fits to 
a grid of synthetic RAVE pipeline metallicities ([m/H] rave) have uncertainty of 

~ 0.2 dex, derived from an external calibration of high-quality spectra of predominantly 
metal- rich stars, which dominate the sample (Zwitter et al. 2008). The grid of synthetic 
spectra utilized extends to [m/H] = —2.5, requiring additional analysis below this limit: we 
provide that here. 



2 The quantity [m/H] here designates the total amount of metals with respect to the solar value. We will 
use [Fe/H] to refer to only the iron content. 
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2.2. Isolating candidate Very Metal Poor Stars in RAVE 

The RAVE survey is designed primarily to obtain excellent line-of-sight velocities. This 
has been achieved, with accuracy and precision of a few km/s for the vast majority of the 
stars (Steinmetz et al. 2006). A large subset of the spectra has sufficiently high signal- 
to-noise ratio for the RAVE pipeline also to provide good quality estimates of the values 
of stellar parameters (Zwitter et al. 2008). At the start of the present investigation, the 
RAVE (internal) database contained over 147,000 spectra from which stellar parameters 
were derived. How do we identify VMP stars among these? Stars with true metallicity near 
or below the lower limit of the synthetic spectra in the RAVE pipeline (—2.5 dex) will be 
assigned metallicities close to that lower limit. A metallicity cut in the RAVE database of, 
say, 3cr above the lower bound, [ul/H]rave< —1.8, should then isolate probable VMP stars. 
Adopting this cut and isolating only normal stars gives a sample of 622 potential metal-poor 
stars. 

Furthermore, extremely metal-poor stars will be erroneously assigned to significantly 
higher metallicities by the RAVE pipeline. The RAVE spectrum of a typical VMP star 
is nearly featureless, with only the Ca Triplet lines being measurable, as shown in Fig- 
ure [TJ There exists a degeneracy between the stellar parameters and metallicity, in that the 
strengths of the Ca Triplet lines can be equal for both a lower-metallicity cooler giant and 
a higher-met allicity hotter turn-off star. A temperature difference of 100 K corresponds to 
a [m/H] rave difference of about 0.1 dex, assuming the parameters of the star still lie along 
an appropriately metal-poor old evolutionary track. Hence stars for which the RAVE T c g 
value is significantly hotter than the value that is derived by other indicators are likely to be 
stars whose real [m/H] values are significantly lower than [hi/H]rave- This expectation is 
confirmed by identification in the RAVE database of several known VM P stars, incorrectly 



allocated [m/H] rave > -1.8, including CD -38°245, with [Fe /H] = -4.2 flCavrel et al.ll2004l ) 



and CS 29502-092 with [Fe/H] = - 2.76 flAoki et all 120071 ). As anticipated, known VMP 



stars with anomalously high [m/H] rave values also have anomalously high estimated stellar 
effective temperatures from the RAVE pipeline. 

Temperature anomali es may be identifie d using photometric temperatures. Photometry 



from the 2MASS survey (ICutri et al.l 120031 ) is available for the RAVE sample, so that an 



effecti ve temperatur e , T p h n t, may be calculated using the (de- reddened) color-T phot relation- 



ship of lAlonso et al.l ( 119991 ). with due care to use the 2-color relations to check for cool binary 
companions. The mean value of E(B-V) for our candidate metal-poor stars is 0.11 mag (me- 
dian 0.06 mag). We assume conservatively that the stars of interest are beyond the reddening 
layer, since their estimated distances based on the stellar parameters are typically several 
hundred parsecs to a few kiloparsecs from the Sun. 
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These T p hot values allow us to adopt the following modified criterion to define a candidate 
VMP star: 



[m/H] RAVE - 0.1((T RAVE - T phot )/(100K)) < -1.8. 

Application of this selection criterion to the RAVE database of non-variable stars within 
the appropriate temperature range to be halo turnoff/red giant stars (7000 < T R ave < 
3800 K) yields 320 stars in addition to those selected by the straight cut [m/H] R AVE< —1.8. 
The total number of stars we selected for further analysis of their RAVE spectrum, based 
on the above criteria, is then 942. 



3. Derivation of [Fe/H] Abundances 

3.1. Equivalent- Width Analysis of RAVE Spectra 

For each candidate VMP star we measure the equivalent width (EW) of the Ca Triplet 
and Fe I lines from the RAVE spectrumjfl These EW va lues are then the input for an 
abundance analysis using the MOO G program (jSneden 1973 ). assuming L TE and 1-D, plane 



parall el Kurucz model atmospheres^ We adopted the NLTE corrections of Starkenburg et al. 



(120101 ) for all our stars (in the end significant only for stars with [Fe/H] < —2.5). 



The estimates of stellar gravity from the RAVE pipeline analysis identify most VMP 
candidates as giants: 86% of these stars have RAVE log g < 3, with the mean RAVE log g 
value being 1.0 ± 0.9. Stellar gravity values from the RAVE pipeline are accurate to better 
than 0.5 dex for the typical metal-rich stars in the RAVE database (Zwitter et al. 2008). 
For metal-poor stars the lower metallicity cut-off of the synthetic spectrum grid biases the 
pipeline analysis to overestimate, systematically, both T eff and log g values, so the actual 
fraction of true giants in the VMP sample should be even higher. Indeed, there were only 
four stars whose location on a reduced proper-motion diagram indicated they were dwarfs. 
For these we adopt log g = 4.5. The stellar surface gravit y for the remaining st ars was 



obtained by fitting to an old (12 Gyr) Yonsei-Yale isochrone ( iDemarque et al.l 120041 ) with a 
metallicity matching the current metallicity estimate of the star (the procedure is iterated). 
This gravity value was adopted in the calculation of the model atmosphere (while the Yonsi- 



3 At least one of the Ca II lines is too strong in the majority of the candidates. All three Ca II lines could 
be used in 132 stars; for these the mean dispersion is 0.26 dex, consistent with our final estimated error. 



4 The most recent versions of the Kurucz atmospheres can be found at http://kurucz.harvard.edu/ 
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Yale isochrones do not extend beyond the RGB, < 0.1 dex error in final abundance results 
should a star be actually on the Asymptotic Giant Branch). Our analysis was iterated until 
the derived (from CaT) metallicity did not differ by more than 0.1 dex from the value used in 
calculating the model atmosphere and isochrone fit. We were able to derive Ca abundances 
for 771 stars out of our 942 star initial selection from RAVE: these 771 stars form our VMP 
candidate sample. 

The next step is to provide an external calibration of our abundance scale, derived from 
high signal-to-noise echelle spectra. 

3.2. Metallicity Calibration via Echelle Observations 

We obtained echelle spectroscopic data for 112 candidate metal-poor stars selected 
from the RAVE database using a preliminary version of the criteria described above. The 
observations were carried out with several telescope/spectrograph combinations, includ- 
ing Magellan-Clay/MIKE, APO-3.5m/ARCES, AAT/UCLES, Max Planck-2.2m/FEROS, 
CFHT/Espandons and VLT/UVES. Full details of the data aquisition and reduction are de- 
ferred to our paper on the derived elemental abundances (Fulbright et al, in prep). Briefly, 
all spectrographs delivered a resolving power greater than 30,000 and, with the exception 
of the UCLES set-up, the wavelength region covered from below 4000 A to beyond 8000 A, 
albeit with some coverage gaps. The effective wavelength range for the UCLES spectra is 
4460-7260 A. In each case the data were reduced using standard methods for echelle data, 
utilizing pipeline programs when available. The typical S/N level of the spectra is greater 
than 100/pixel and often exceeds 200/pixel. 

The abundance analysis utilized Kurucz stellar atmospheres and the MOOG program, 
with now the metallicity of the stellar atmosphere for each star set equal to the value of 
[Fe II/H] derived from the analysis. The value of the stellar effective temperature was 
derived using the excitation temperature method based on Fe I lines. The mean difference of 
T cx — T phot is —66 ± 138 K. The value of the stellar gravity was again taken to be that value 
corresponding to the derived T e ff on an old (12 Gyr) Yonsi-Yale isochrone of the appropriate 
metallicity. The microturbulent velocity was set to the value that minimized the slope of the 
relationship between the iron abundance derived from Fe I lines and the value of the reduced 
equivalent-width of each line. The Fe I lines measured in the RAVE spectra are usually too 
weak in true VMP stars to give reliable results. The iron abundances derived from these 
lines were useful, however, in rejecting from the analysis those (few) stars for which the Ca 
Triplet lines appeared weak, but the Fe I lines were strong. Further inspection of the spectra 
of these stars revealed that the Ca lines are likely weakened by emission cores. 
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The iron abundances from the echelle data provide an immediate check on our selection 
criteria. The result is encouraging: of the 92 stars with echelle data for which we predicted 
[Fe/H] < -2 from the RAVE data, 87 (95%) indeed have [Fe/H] values below -2. 

Calibration of the RAVE equivalent width abundances is then achieved through least- 
squares fitting between the echelle-based iron abundance, [Fe/Hjni-Res, and the equivalent- 
width based calcium abundance from the RAVE spectra, [Ca/H] R AVE- The correlation be- 
tween these two quantities is shown in Figure H] and the relationship is: 

[Fe/H] ffi _Rcs = 0.93[Ca/H] RAVE - 1.33, 
with a correlation coefficient of 0.82, and standard deviation of the residuals of 0.25 dex. 



4. The Stellar [Fe/H] Abundance Distribution Function at low metallicity 

The calibration obtained from the echelle data was applied to the 771 VMP candi- 
dates from RAVE for which the MOOG analysis provided the calcium abundance. The 
RAVE-based [Fe/H] abundance distribution function that resulted is shown as the solid-line 
histogram in Figure [3J this contains a total of 612 stars with true iron abundance [Fe/H] 
< -2 dex. 



4.1. Distribution Function Completeness 



Given that our candidate sample had an upper metallicity cutoff determined by uncer- 
tain quantities, significant incompleteness at the high-metallicity end is expected. We check 
this using published catalogs of candidate VMP stars and a second RAVE study. 

A SIMBAD-aided literature search showed that only 47 of our 612 very metal-poor stars 
([Fe/H] < —2 dex) have previously been pr oposed as be i ng me tal-poor, from low-re s olutio n 
spectroscopic analyses, including 21 from iFrebel et al.l (120061 ) and 11 from iBondl (|1980l ). 
Ten of the 138 stars for which we derive —2 < [Fe/H] < — 1 also have entries in SIMBAD 
identifying them as VMP candidates. 

The most uniform compa rison is with the Hamb urg-ESO survey (HES) catalog of over 
20,000 candidate VMP stars Jchristlieb et al.lboosh . Cross-identification with the RAVE 
database (at that time with 200,000 entries) yielded 473 matches, with HES metallicity 
estimates for 296 stars, all with [Fe/H]nES < —2. 109 of these stars are included in our 
RAVE VMP sample. The RAVE pipeline gave [ul/H]rave< — 1 for 76 of the remaining 187 
stars. We reanalyzed these RAVE spectra as above, obtaining [Fe/H] < — 1 for only 22 stars, 
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with just 4 with [Fe/H] < —2, of which the most metal-poor has [Fe/H] = —2.26. These 
stars are therefore all very close to the VMP threshold. 



The high-resolution sample of iRuchti et al.l ( 120101 ) was selected from the same RAVE 



catalog as this study and includes 20 stars for which those authors derive [Fe/H] < —2 
from their echelle spectra but which were not selected by our VMP criteria. We re-analyzed 
the RAVE spectra of these stars following the procedures developed above and found good 
agreement: [Fe/H] < —2.0 for 15 stars, and all 20 stars having [Fe/H] < —1.7. The lowest 
iron abundance of this group is —2.41 dex, within 2a of our calibration. Approximately half 
of these stars have [m/H] rave between —1.5 and —1.8, the remainder having RAVE spectra 
of low quality. As expected, measuring errors in our parent RAVE sample lead to increasing 
incompleteness in our final sample, at above [Fe/H]~ —2.5. 



4.2. The low metallicity tail of the distribution function 

At very low metallicities, we test for true VMP stars observed by RAVE, but excluded 
from our sample, by cross-matching the RAVE database with publi shed catalogs of con - 



firmed VMP stars. There are 253 RAVE stars in the HERES project (IBarklem et al.ll2005l ). 
which provides high-resolution follow-up for HES catalog stars. There are only two candidate 
VMP matches: J112243.4-020936 (RAVE) is HE1120-0153 (HERES), and J132244.1-135531 
(RAVE) is HE1320-1339 (HERES). Our VMP star criteria identified both, and our calibra- 
tion gave each an iron abundance of [Fe/H] = —2.6. The HERES values are [Fe/H] = —2.77 
and —2.78, respectively, in excellent agreement. No confirmed very metal-poor star has been 
excluded from our RAVE sample. 

There are three stars in our RAVE sample at or below —4 dex. These are C0022448- 
172429 for which we derive, from this calibration, [Fe/H] = —4.0 dex, CD — 38°245 for 
which we derive [Fe/H] = —4.2 dex, and a third star, with a RAVE spectrum at our min- 
imum acceptable RAVE S/N cutoff, which our calibration gives [Fe/H] = —4.0 dex. Our 
echelle-based result for C0022448-172429 is [Fe/H] Hi -R C s = -4.02, and for CD -38°245 is 
[Fe/H] Hi-Rcs = —4 .20. CD — 3 8 °245 was previously confirmed from echelle data to have 



[Fe/H] = —4.2 by ICayrel et al.l (120041 ) . Our RAVE-based calibration is indeed providing 



true iron abundances to ~ 0.25 dex accuracy, even at the lowest known metallicities. 

Our analysis is the first confirmation of the 'ultra metal-poor' nature, [Fe/H] < —4 dex, 
of C0022448-172429. The low resolution, objective prism HES spectrum for C0022448-172429 
(Christlieb et al. 2008, where this star is identified as HE0020-1741) provided estimates of 
—3.0 and —3.3, depending on which method those authors used. Christlieb et al. further 
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estimate that C0022448-172429 is carbon-rich, with [C/Fe] = +1.0. Our echelle analysis 
confirms the carbon-rich nature, with our value being [C/Fe] = +0.9, and we find this star 
has an extremely high nitrogen abundance, [N/Fe] = + 3.2. Thus this star joins the select 



group of carbon-enhanced metal poor (CEMP) stars (see lNorris. Gilmore. Wyse. et aUl2010 
for a recent discussion of this class of object). 

Although this makes our selection function not internally consistent, adding the iden- 
tified extra VMP stars after re-analysis of their RAVE spectra (15 from Ruchti et al., plus 
4 from the HES cross-check) to our main sample of 612 VMP stars, results in the dashed 
histogram in Fig. [3j This contains 631 stars with [Fe/H] < —2 dex, i.e 'Very Metal- Poor', 
and two stars with [Fe/H]< —4, i.e. 'Ultra Metal-Poor' and one at ~ —4 dex. 567 of these 
VMP stars were not known to be VMP previously, and one of the two UMP stars is a new 
confirmation. The remaining star at ~ —4 dex has no previous published abundance deter- 
mination and due to its marginal quality RAVE spectrum remains as a candidate, awaiting 
scheduled high resolution follow-up; if confirmed this would be a new discovery. 



5. Conclusions 

The RAVE Survey allows identification of large and complete samples of apparently 
bright, very metal poor (VMP) stars, ideal for subsequent detailed analysis. We have defined 
a sample of 612 stars with [Fe/H] < —2, based on some 200,000 spectra in the RAVE database 
(which currently contains 350,000 spectra). This includes two stars with [Fe/H] < —4 dex, 
one of which was previously known, the other of which is a new ultra metal-poor star, with 
the previous estimate of iron abundance being too high by 0.7 to 1.0 dex. This star belongs 
to the interesting carbon-enhanced metal poor class, having [C/Fe] = +0.9, and [N /Fe] 



3.2. A third (new) metallicity —4 candidate remains to be confirmed. Only six (IFrebel 



20101 ) such ultra metal-poor stars were known prior to this study. Comparison to other 
samples indicates that our completeness limit at the metal-rich end is [Fe/H] « —2.5. It is 
very likely our sample is complete for more metal-poor stars. 

We note here, and will discuss elsewhere, that the shape of our metallicity distribu- 
tion function below [Fe/H]~ —2.8, where we are complete, is in excel lent agreement with 



the independe nt determinations from the (low-resolution) HES survey (iSchorck et al.l 12009 



Li et al. 



2010|). 



Our technique uses a selection criterion based on the values of [m/H] rave and effective 
temperature from the RAVE pipeline analysis, plus a photometric effective temperature. 
The original RAVE spectrum of a star that passes the criterion is re-analysed to determine 
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the equivalent widths of the Ca T lines and Fe I lines. We calibrate the (re-measured) 
RAVE metallicity onto true iron abundance, [Fe/H], through an echelle-based calibration 
derived in this paper. The efficiency of our technique for finding VMP stars is very high, 
demonstrated by the fact that ~ 95% of the stars that we determine by our calibration 
to have [Fe/H]< —2 indeed have echelle-based [Fe/H] in this range. There is no need for 
follow-up, higher spectral resolution data to confirm this derived iron abundance. Since the 
mean I magnitude of the RAVE sample of VMP stars is ~ 10.5, follow-up detailed elemental 
abundance analyses are both straightforward and efficient. The RAVE survey is on-going, 
and the final RAVE database should at least triple the current sample size. 

Facilities: UKST (6dF), ARC (ARCES), CFHT (Espandons), MagellamClay (MIKE), 
VLT:Kueyen (UVES), Max (FEROS) 
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Fig. 1. — RAVE spectra of three metal-poor stars, with our derived iron abundance estimates 
(consistent with previous estimates in the literature). Note the dominance of the Ca II Triplet 
lines, and the growing strength of Fe I and other species as [Fe/H] increases. 
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Fig. 2. — Calibration of the calcium abundance derived from the RAVE spectra onto iron 
abundance, using [Fe/H] from high-resolution echelle spectra, for 112 stars. 
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Fig. 3. — Distribution of calibrated iron abundances [Fe/H], for the sample of VMP stars 
selected from the RAVE database (solid line). The dashed histogram includes the VMP 
stars from the re-analysis of RAVE spectra for candidates from the HES survey and from 
Ruchti et al. (2010). 



